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Abstract 

The production of -D*^ mesons in deep inelastic ep scattering has been meas- 
ured for exchanged photon virtuahties 5 < < 1000 GeV^, using an integrated 
luminosity of 363 pb~^ with the ZEUS detector at HERA. Differential cross 
sections have been measured and compared to next-to-leading-order QCD cal- 
culations. The cross-sections are used to extract the charm contribution to the 
proton structure functions, expressed in terms of the reduced charm cross sec- 
tion, cr^gjj. Theoretical calculations based on fits to inclusive HERA data are 
compared to the results. 
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1 Introduction 



The measurement of charm production in deep inelastic ep scattering (DIS) is a powerful 
tool to study quantum chromodynamics (QCD) and the proton structure. In leading-order 
QCD, charm production occurs through the boson-gluon fusion (BGF) process '~f*g — cc, 
which is directly sensitive to the gluon content of the proton. Different approaches to 
the calculation of the heavy-quark contribution to the proton structure functions are 
currently used in global analyses of parton density functions (PDFs) [IHl]- Comparisons 
to measurements of charm production in DIS provide direct tests of these approaches [5] . 
It has also been shown recently that a combined analysis of charm production and inclusive 
DIS data can provide a competitive determination of the charm-quark mass [SH?]. 

Several measurements of charm production in DIS have been performed at HERA, exploit- 
ing reconstructed [8], [SHTO] and D*^ [TTl - fTT] mesons, semi-leptonic decays [T8] . 
and inclusive lifetime methods [T9l|20] to tag charm. In this paper, a new high-statistics 
measurement of D*^ production via the reaction 

e{k)p{P) ^ e'{k')D*^{p^')X 

is presented. The symbols in parenthesis represent the four-momenta of the incoming 
(k) and outgoing electron {k'), of the incoming proton (P), and of the produced D*^ 
(p^*). The measurement is performed for photon virtualities, = —q^ = —{k' — k)"^, 
in the range 5 < < 1000 GeV^ and for inelasticities, y = {P ■ q)/{P ■ k), in the range 
0.02 <y < 0.7. 

The D*~^ meson^ were reconstructed through the decay D*^ — )■ D^ti^ with £)" — )■ K~ti^ . 
Differential cross sections are presented as a function of Q^, the Bjorken-x variable, 
and of the fraction of the exchanged-photon energy transferred to the D*'^ meson in the 
proton rest frame, z^* = (P ■ p^*)/{P ■ q), as well as of the D*^ pseudorapidity, r]^* , and 
the transverse momentum, p-^* , in the laboratory frame B 

Double-differential cross sections in and y are presented and used to extract the charm 
contribution to the proton structure functions in the form of the reduced charm cross 
section, cr^g^. Previous measurements and theoretical calculations are compared to the 
results. 



^ Hereafter the charge conjugated states are implied. 

^ The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the 
proton beam direction, referred to as the "forward direction", and the X axis pointing towards the 
centre of HERA. The coordinate origin is at the nominal interaction point. The pseudorapidity is 
defined as 77 = — In (tan|), where the polar angle, 0, is measured with respect to the proton beam 
direction. 
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2 Experimental set-up 



The measurement was based on e p collisions collected with the ZEUS detector at HERA 
in the period 2004-2007 with an electroij^l beam energy, E^, of 27.5 GeV and a proton 
beam energy, Ep, of 920 GeV, corresponding to a centre-of-mass energy -y/i = 318 GeV. 
The corresponding integrated luminosity, C = 363 ±7 pb~^, is four times larger than that 
used for the previous ZEUS measurement [TT] . 

A detailed description of the ZEUS detector can be found elsewhere [2T]. In the kinematic 
range of the analysis, charged particles were tracked in the central tracking detector 
(CTD) [22] and in the microvertex detector (MVD) [23]. These components operated 
in a magnetic field of 1.43 T provided by a thin superconducting solenoid. The CTD 
consisted of 72 cylindrical drift chamber layers, organised in nine superlayers covering the 
polar-angle region 15° < 9 < 164°. The MVD consisted of a barrel (BMVD) and a forward 
(FMVD) section with three cylindrical layers and four vertical planes of single-sided silicon 
strip sensors in the BMVD and FMVD respectively. The BMVD provided polar-angle 
coverage for tracks crossing the three layers from 30° to 150°. The FMVD extended the 
polar- angle coverage in the forward region down to 7°. For CTD-MVD tracks that pass 
through all nine CTD superlayers, the momentum resolution was o-{pt)/pt = 0.0029prffi 
0.0081 © 0.0012/pr, with pr in GeV. 

The high-resolution uranium-scintillator calorimeter (CAL) [2^ consisted of three parts: 
the forward, the barrel, and the rear (RCAL) calorimeters. Under test-beam conditions, 
the CAL single-particle relative energy resolutions were a{E)/E = 0.18/a/E for electrons 
and a{E)/E = 0.35/^/E for hadrons, with E in GeV. The energy of electrons hitting 
the RCAL was corrected for the presence of dead material using the rear presampler 
detector [25] and the small angle rear tracking detector (SRTD) [26] . 

The luminosity was measured using the Bethe-Heitler reaction ep — ?■ e'jp by a luminosity 
detector which consisted of two independent systems: a lead-scintillator calorimeter [27] 
and a magnetic spectrometer [28] . 

3 QCD calculations 

Cross sections for heavy-quark production in DIS were calculated at next-to-leading or- 
der (NLO), i.e. O(a^), in the fixed-flavour-number scheme (FFNS), in which only light 
flavours and gluons are present as partons in the proton and heavy quarks are produced 

^ Hereafter "electron" refers to both electrons and positrons unless otherwise stated. 



2 



in the hard interaction [29j. The program HVQDIS [301 [3T] was used to compute single- 
and double-differential D*~^ cross sections. 

The parameters used as input to HvQDis are listed below, together with the variations 
used to evaluate the uncertainty on the theoretical prediction: 

• charm-quark pole mass: rric = 1.50 ± 0.15 GeV; 

• renormalisation (/i/j) and factorisation {^f) scales: = = + 4m^, varied 
independently up and down by a factor two; 

• strong coupling constant in the three-flavour FFNS: af^^{Mz) = 0.105 ± 0.002; 

• the PDFs and their uncertainties, taken from a FFNS variant [3] of the HERAPDFl.O 
fit [32]. The central PDF set was obtained from a fit performed using the same values 
of rric, /i_R, fiF and as used in the HvQDis program. For each variation of these 
parameters in HvQDis, a different PDF set was used, in which the parameters were 
varied consistently. 

The NLO calculation provided differential cross sections for charm quarks. The fragment- 
ation model described in a previous publication [5] was used to compare to the measured 
D*^ cross sections. This model is based on the fragmentation function of Kartvelishvili et 
al. [33], controlled by the parameter ax, to describe the fraction of the charm momentum 
transferred to the D*^ mesons. It also implements a transverse fragmentation component 
by assigning to the D*~^ meson a transverse momentum, kr, with respect to the charm- 
quark direction. The uncertainty on the fragmentation model was estimated by varying 
aK and the average kx according to the original prescription [5]. The fraction of charm 
quarks hadronising into D*~^ mesons was set to /(c — ?■ D*^) = 0.2287 ± 0.0056 [M] . 

For the inclusive cross section, theoretical predictions were also obtained in the generalised- 
mass variable- flavour- number scheme (GM-VFNS). In this scheme, charm quarks are 
treated as massive particles for < ml and as massless partons for 3> ml, in- 
terpolating in the intermediate region [53HTF] . The calculation was performed using the 
Roberts-Thorne (RT) "standard" [38|l39] variant of the GM-VFNS at NLO, correspond- 
ing to 0{al) for the < ml part and to 0{as) for the ^ ml part. PDFs obtained 
from the HERAPDF1.5 [IQ] fit to inclusive HERA data were used. The central prediction 
was obtained for = 1.5 GeV. To evaluate the theoretical uncertainty, the calculation 
was repeated varying the PDF set and its parameters according to the systematic vari- 
ations associated with the HERAPDF1.5 fit. The dominant source of uncertainty was the 
charm-quark mass, which was varied in the range 1.35 < m^ < 1.65 GeV. 
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4 Monte Carlo samples 



Monte Carlo (MC) samples were used to calculate the experimental acceptance and to 
estimate the background contamination. MC samples of charm and beauty DIS events 
were generated using Rapgap 3.00 [H]. The main sample consisted of events generated 
according to the LO BGF process. Radiative QED corrections to the BGF process were 
included through Heracles 4.6 |42]. Additional Rapgap samples were generated for 
diffractive charm production and for the resolved-photon processes gg — )■ cc and eg — )■ eg, 
in which one of the incoming partons originates from the exchanged photon. Charm 
photoproduction was simulated using Pythia 6.2 [13]. 

Both Rapgap and Pythia use parton showers to simulate higher-order QCD effects and 
use the Pythia/ Jetset hadronisation model [13]. All samples were generated using 
the CTEQ5L [H] proton PDFs and, for resolved-photon processes, the GRV-G LO [15] 
photon PDFs. The diffractive samples were generated using the "HI fit 2" jlB] diffractive 
PDFs. The heavy-quark masses were set to ttIc = 1.5 GeV and = 4.75 GeV. Masses, 
widths and lifetimes of charmed mesons were taken from PDG2010 |47j . 

The MC samples correspond to about four times the luminosity of the data and were 
passed through a full simulation of the ZEUS detector based on Geant 3.21 [IS]. They 
were then subjected to the same trigger criteria and reconstructed with the same programs 
as used for the data. 

5 Event selection and signal extraction 
5.1 DIS event selection 

A three-level trigger system was used to select DIS events online [2Tlll9l[50] by requiring 
electromagnetic energy deposits in the GAL at the first level and applying loose DIS 
selection criteria at the second and third levels. 

Offline, the hadronic system was reconstructed using energy-flow objects (EFOs) [51] 
which combine tracking and calorimeter information. The electron was identifled us- 
ing a neural-network algorithm [52]. The kinematical variables Q^, y, and x were re- 
constructed using the S method [53]. The variable z^* was reconstructed according to 
z^* = {E^* —p^'')/{2EeyjB), where yjB is the inelasticity reconstructed with the Jacquet- 
Blondel method [51] and E^* and p^* are the D*^ energy and longitudinal momentum, 
respectively. 

The following criteria were applied to select DIS events [53] : 
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• Ee' > 10 GeV, where Egi is the energy of the scattered electron; 

• Uf. < 0.7, yjB > 0.02, where ?/e is the inelasticity reconstructed from the scattered 
electron; 

• 40 <£' — < 70 GeV, where E — Pz is the global difference of energy and longitudinal 
momentum, obtained by summing the electron and the hadronic final state, which is 
expected to be 2Ee = 55 GeV for fully contained events; 

• the Z position of the primary vertex, Z^t^, was required to be in the range \Z^t^\ < 
30 cm; 

• the impact point of the scattered electron on the RCAL was required to lie outside a 
square region around the beam-pipe hole: \Xe\ > 15 cm or |Ye| > 15 cm; 

• 5 < < 1000 GeV^, where is reconstructed with the S method. 

5.2 Selection of D*+ candidates and signal extraction 

The D*~^ mesons were identified using the decay channel D*^ — t- D^tt^ with the sub- 
sequent decay — K~7r~^, where tt^ refers to a low-momentum ("slow") pion accom- 
panying the D^. 

Tracks from the D*~^ decay products were required to have at least one hit in the MVD 
or in the inner superlayer of the CTD and to reach at least the third superlayer. Tracks 
with opposite charge and with transverse momentum p!^'^ > 0.4 GeV were combined in 
pairs to form D° candidates. The track parameters were improved by fitting the two 
tracks to a common vertex. Pairs incompatible with coming from the same decay were 
removed by requiring a distance of closest approach of the two tracks of less than 1 mm, 
and the of two-track vertex fit smaller than 20 for one degree of freedom. The 
tracks were alternately assigned the kaon and pion mass and the invariant mass of the 
pair, M{Ktt), was calculated. Each additional track, with charge opposite to that of the 
kaon track and a transverse momentum > 0.12 GeV, was assigned the pion mass and 
combined with the -D° candidate to form a D*+ candidate. The iTg track was then fitted 
to the primary vertex of the event, obtained exploiting the other tracks reconstructed in 
the event and the constraint from the average position of the interaction point [8]. The 
mass difference AM = MIKtt-Ks) — M{Ktt) was used to extract the D*~^ signal. The D*~^ 
candidates were required to have 1.80 < M{Ktx) < 1.92 GeV, 143.2 < AM < 147.7MeV, 
1.5 <p^* < 20 GeV and < 1.5. 

The distribution of M{Ktt) for D*~^ candidates, without the requirement on M{K7r), is 
shown in Fig. [H Also shown is the distribution of wrong-sign (WS) candidates, obtained 
by combining two tracks with the same charge. The WS distribution provides an estimate 
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of combinatorial backgrounds. A clear peak at the mass is visible in the correct-sign 
(CS) distribution. The excess of CS candidates at masses below the D° peak is due to 
partly-reconstructed decays, mostly — )■ K^n^ii^. 

The distribution of AM for D*^ candidates, without the requirement on AM, is shown 
in Fig. [21 A clear D*~^ peak is seen. The D*~^ signal was extracted by subtracting the 
background estimate from the number of candidates in the signal window 143.2 < AM < 
147.7 MeV. The background estimate was obtained by fitting simultaneously the CS and 
WS distributions to the parametrisation 



where A, B, C, D are free parameters of the fit [56] and C, = AM — m^+. The fit was 
performed in the region AM < 168 MeV. The region with a possible signal contribution, 
140 < AM < 150 MeV, was removed from the fit to the CS distribution. The parameter 
D, which represents the normalisation of the CS background with respect to the WS 
distribution, is slightly larger than unity, D = 1.021 ± 0.005. This is consistent with 
the MC estimation of the additional combinatorial background component in the CS 
distribution due to real — Kn decays associated with a random track to form a CS 
D*+ candidate. The total signal is N^^^!^ = 12893 ± 185. 

The amount of signal lost due to the tails of the mass peak leaking outside the M{Ktt) 
window was estimated by enlarging the mass window to 1.7 < M{K'k) < 2.0 GeV. The 
fraction of additional D*~^ found within the enlarged window was 13%, including the 
contribution from partly reconstructed D^. This fraction, as well as its dependence onp^* 
and rj^* and on the width of the M[Kti) window, was found to be well reproduced by MC. 
The signal in the tails of the D*^ peak outside the AM window was estimated similarly, 
enlarging the signal window to 140 < AM < 150 MeV. The fraction of additional D*^ 
was 6% on average, with a dependence on the transverse momentum of the slow pion, 
due to the momentum and angular resolution degrading at low . This effect is not 
completely reproduced by the MC. An acceptance correction [SS] dependent on pj" was 
then applied, ranging from ^ 10% at = 0.12 GeV to ~ 1% at large p^" . 



6 Cross-section extraction 

The differential cross sections, da^is/dC,, for producing a D*'^ in the "visible" phase space 
1.5 < p^* < 20 GeV, \r]^*\ < 1.5, 5 < < lOOO GeV^ and 0.02 < y < 0.7 was obtained 



WS: ^s(C) 

CS: /es(C) 



^/ws(C), 



as 




■a 



- A- BR -C 
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where A'^^.^^ is the signal extracted in a bin of a given variable ^, A'^ is the pho- 
toproduction background, A,^ is the bin size, A is the acceptance, BR = B{D*^ — )• 
D°7r+) X B{D^ A'~7r+) = 0.0263 ± 0.0004 [57] is the branching ratio, £ is the integ- 
rated luminosity and Cr is the QED radiative correction. 

The background from charm photoproduction {Q^ < 1.5 GeV^) was evaluated using 
the photoproduction MC sample, normalised to the luminosity using the cross sections 
previously measured by ZEUS [SB]. 

The acceptance. A, was calculated as the ratio between the number of reconstructed 
and generated D*^ in the bin, using a signal MC based on a mix of charm and beauty 
production. The beauty MC was normalised to 1.6 times the cross section given by 
Rapgap, consistent with ZEUS measurements [TS IISUH^ . The charm MC contained non- 
diffractive and diffractive components, summed according to the relative cross sections as 
given by Rapgap. The normalisation of the charm MC was adjusted such that the sum 
of all the MC components reproduced the number of D*~^ mesons in the data. Resolved- 
photon processes were not included. They were only used for systematic checks. The rj^* 
and p^* distributions of the charm MC were reweighted [SSj to improve the agreement 
with data, with the p-^* weights dependent on Q^. 

The acceptance as determined by the MC was corrected to account for imperfections in 
the simulation of the trigger and track-reconstruction efficiencies. One of the main sources 
of track-reconstruction inefficiency for charged pions and kaons were hadronic interactions 
in the material between the interaction point and the CTD. This effect was studied using 
special tracks from ep — )■ ep^ with p° — vr+vr" events, reconstructed from MVD hit 
information alone [52] • For these tracks, an extension into the CTD was searched for. 
In addition, the pt dependence of the tracking efficiency was studied by exploiting the 
isotropic angular distribution of pions from Kg decays. The studies showed that the MC 
slightly underestimated the effect of nuclear interactions. For central pions with pj- ^ 
1 GeV, the track-reconstruction inefficiency due to hadronic interactions was measured to 
be (7± 1)% while the MC predicted 5%. The track-efficiency correction was applied as a 
function of rj and pt of each track. For px > 1-5 GeV, no correction was necessary. 

The acceptance ranges from A ~ 10% in the lowest p!^ and bins to ^ ^ 45% in the 
highest p^* and bins. Fig. [3] shows N^*^^/A^ for = pfj^' , r]^' ,Q^, y and z^* . The 
sum of the different MC samples is compared to the data. The agreement is satisfactory. 

The cross sections were corrected to the QED Born level, using a running coupling con- 
stant aem(Q^), such that they can be compared directly to the QCD predictions from 
the HvQDis program. The radiative corrections were obtained as Cr = o'vi^'^/'^vtif ' where 
cr^is'''^ is the Rapgap cross section with the QED corrections turned off but keeping ctem 
running and cr^^jf is the Rapgap cross section with the full QED corrections, as in the 
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standard MC samples. 



7 Systematic uncertainties 

The experimental systematic uncertainties are listed below [SS], with their typical effect 
on the measured cross sections is given in parenthesis: 

6i energy-scale uncertainty on the hadronic system of ±2% (±1%, up to ±10% at low 

y); 

52 electron energy-scale uncertainty of ±1% [63] (±1%, up to ±7% at low y); 

53 alignment uncertainty on the electron impact point on the RCAL, estimated by varying 
the cut on the electron position in the MC by ±2 mm separately for the Xe and Ye 
coordinates [M] (±7% at low and low y, negligible at large Q^); 

54 uncertainty on the position of the electron impact point on the RCAL due to imperfec- 
tions in the simulation of the shower shape and of the detector resolution, estimated 
by loosening the cut on the electron position by 1 cm {\Xe\ > 14 cm or \Ye\ > 14 cm) 
both in data and in MC (up to ±10% at low y and low Q"^, negligible at large Q^); 

^5 uncertainty on the background shape in AM, estimated by replacing the function 
/cs(C) by /4(C) = AC§ +BC + CC^ (±0.3%); 

^6 a further uncertainty on the background shape, evaluated by reducing the fit range 
from AM < 168 MeV to AM < 165 MeV (±0.5%); 

^7 uncertainty on the amount of signal outside the AM window, evaluated by varying 
the pj'^- dependent correction by its uncertainty (±1.5%, up to ±3% at low pr); 

^8 uncertainty on the amount of signal outside the M{Kt[) window, estimated by com- 
paring data and MC in an enlarged mass range (±2%); 

^9 uncertainty on the track-efficiency, evaluated by varying the track efficiency correction 
applied to MC by the associated uncertainty (±2%); 

5iQ uncertainty on the trigger efficiency, evaluated using independent triggers (±0.5%); 

511 statistical uncertainty on the calculation of the acceptance (±1%); 

512 uncertainty on the normalisation of the beauty MC of ±50% to cover the range allowed 
by ZEUS measurements [59l[6T] (±0.3%); 

513 uncertainty on the normalisation of the photoproduction MC of ±100% (up to ±3% 
at high y, but negligible elsewhere); 
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5i4 uncertainty on the normalisation of the diffractive charm MC of ±50% to cover the 
range allowed by data-MC comparison and by previous ZEUS results [51] (up to 
±4.5% at low y, but negligible elsewhere); 

(5i5 uncertainty due to the resolved-photon component, evaluated by adding the resolved- 
photon samples to the charm MC normalised according to the generator cross section 
(±2%); 

^16 uncertainty on the MC reweighting as a function of p^* and Q^, which was varied by 
±50% (±2%); 

(5i7 uncertainty on the MC reweighting as a function of r]^* which was replaced by a MC 
reweighting as a function of y (from —2% to ±3%, depending on y)\ 

(5i8 uncertainty on the integrated luminosity of ±1.9%; 

(5i9 uncertainty on the branching ratio BR of ±1.5%. 

All the systematic uncertainties, except the overall normalisations ^ig and ^ig, were added 
in quadrature to the statistical uncertainties to obtain the total error bars in the figures. 



8 Results 

Single- and double-differential cross sections have been measured in the phase space 

5 < Q2 < 1000 GeV^; 0.02 <y< 0.7; 1.5 < p^' < 20 GeV; |r/^* | < 1.5. 

Differential cross sections in pi^* , r]^ and z^* are reported in Tables [TH3] and in Fig. IH 
The cross section decreases steeply with pJ^ and is almost constant in rj^ . The NLO 
calculations based on HvQDis and the Rapgap MC implementing the leading-order BGF 
process are compared to the data. As the Rapgap MC is based on leading-order matrix 
elements, it is not expected to estimate the normalisation correctly. Therefore the Rap- 
gap prediction was normalised to the data, scaling it by 1.1, to allow a direct comparison 
of the shapes. The data are well described by the NLO calculation and by Rapgap with 
the exception of the shape in z^* , which is not well reproduced by the NLO calculation, 
suggesting possible imperfections in the fragmentation model. 

Differential cross sections m. Q^, y and x are reported in Tables HHB] and in Fig. O The 
results are reasonably well described by the NLO calculation. The MC predictions repro- 
duce the shapes of the data, except for the high-Q^ tail, where the MC prediction is too 
high, and for da/dy, where the prediction is too low at low y and too high at large y. 
These imperfections in the MC are to be expected in the absence of higher-order terms 
in Rapgap. 
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Visible cross sections in two-dimensional bins of and y, Uvis, are given in Table [71 
The corresponding bin-averaged double-differential cross sections are shown in Figs. [HI 
and [71 The values of the individual systematic uncertainties on the double-differential 
cross sections are given in Table [3 Measurements performed in the same phase space 
by the HI Collaboration [HI [17], which are the most precise previous measurement of 
£)*+ production in DIS, are compared to the present results. The two data sets are 
in agreement and have similar precision. The double-differential cross sections are well 
described by the NLO calculation. 

In a previous ZEUS measurement [TT], a possible excess in the D*"'" yield in e~p collisions 
was observed with respect to e^j) collisions. The ratio of observed rates, increasing with 
was r^'vjr^-^v = 1.67 ± 0.21(stat.) for 40 < < lOOO GeV^ The measurement 
was based on a luminosity of 17 (65) pb^^ of e~p {e~^p) collisions. The present meas- 
urement is based on an independent data set, consisting of 187 (174) pb~^ of e~p (e'^p) 
collisions. Fig. [S] shows the cross-section ratio as a function of Q^. Only statistical uncer- 
tainties are shown since systematic effects mostly cancel in the ratio. No deviation from 
unity is observed, confirming the original interpretation of the e~p excess as a statistical 
fluctuation. 



9 Charm reduced cross sections 

The reduced cross section for charm, a^^^j, and the charm contribution to the proton 
structure functions, ^2^*^ and F^'^, are defined as: 



,2 



<d(x, Q\ s) = Ff (x, Q2) - ^Ffix, g2), 
where F+ = 1 + (1 — y)"^- 

The HvQDis program was used to extrapolate the measured visible D*~^ cross sections in 
bins of y and Q^, (Xvis, to the full phase space: 

_ f beauty^ / " red, HvQDiS l"^' 

''red^'" 



/ ^■z\ ( beauty\ / ^red, HvQDiS Q 



''"vis, HvQDis 



where <J^^^^^ is the beauty contribution as predicted by the Rapgap MC, normalised as 
discussed in Section [SI and a'^l^ hvqdis> '^vIs.hvqdis are the charm reduced and the visible 
D*^ cross sections, respectively, as given by HvQDis. The reference values of x and 
were chosen close to the average x and of the bins. The kinematic acceptance of the 
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visible phase space, defined as Aps = a^is/icr'^'^ ■ 2 f{c — t- D*^)), where a'^'^ is the charm 
production total cross section in the y and bin, ranges from 17% to 64%, depending 
on the bin. 

Following the method used in the previously published combination of ZEUS and HI 
results [5], the HvQDis and fragmentation variations described in Section [3] were used to 
determine the theoretical uncertainty on the extraction of a^^^. The scales fiR and fip 
were varied simultaneously rather than independently as in the theoretical uncertainty 
for the differential cross sections. An additional uncertainty originates from the subtrac- 
ted beauty component that was varied by ±50%. The theoretical uncertainties due the 
extrapolation on cr^g^{x, Q"^) are given in Table O The experimental part of the uncertain- 
ties on Q"^) is defined as the quadratic sum of the statistical and the experimental 
systematic uncertainties described in Section [3 

The results are reported in Table [10] and are shown in Fig. [91 The combined result based 
on previous HI and ZEUS charm measurements [S] and a recent ZEUS measurement 
with mesons [9], not included in the combined results, are also shown. All three 
measurements are in good agreement. The D* measurement has a precision close to that 
of the combined result in some parts of the phase space. The GM-VFNS calculation, 
based on the HERAPDF1.5 parton-density fit to inclusive HERA data, is compared to 
the present measurement and shown in Fig. [TOl The uncertainty on the prediction is 
dominated by the charm-quark mass. The prediction is in good agreement with the data. 



10 Conclusions 

Differential cross sections for the production of D*"^ mesons in DIS have been measured 
with the ZEUS detector in the kinematic range 

5 < < 1000 GeV^; 0.02 <y< 0.7; 1.5 < p^* < 20 GeV; |r/^* | < 1.5, 

using data from an integrated luminosity of 363 pb~^. The new data represents one of the 
most precise measurements of charm production in DIS obtained to date. The data are 
reasonably well described by NLO QCD calculations and are in agreement with previously 
published results. 

The measurements have been used to extract the reduced cross sections for charm cr^g^. 
A GM-VFNS calculation based on a PDF fit to inclusive DIS HERA data agrees well 
with the results. This demonstrates a consistent description of charm and inclusive data 
within the NLO QCD framework. 
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n* 
Pt 

(GeV) 


da 
dp!j^* 

(nb/ GeV) 


fstat 
(%) 


(%) 




1.50 : 1.88 


2.16 


9.9 


+7.0 
-5.5 


1.03 


1.88 : 2.28 


2.30 


5.8 


+5.4 
-5.8 


1.04 


2.28 : 2.68 


1.95 


4.4 


+5.0 
-4.4 


1.03 


2.68 : 3.08 


1.63 


4.0 


+4.7 
-4.0 


1.03 


3.08 : 3.50 


1.22 


3.8 


+4.9 
-4.2 


1.04 


3.50 : 4.00 


9.71x10-1 


3.4 


+4.4 
-3.7 


1.03 


4.00 : 4.75 


6.26x10^1 


3.2 


+4.2 
-3.5 


1.05 


4.75 : 6.00 


3.32x10^1 


3.0 


+4.3 
-3.7 


1.01 


6.00 : 8.00 


1.21x10-1 


4.1 


+4.1 
-3.8 


1.06 


8.00 : 11.00 


3.31x10-2 


6.0 


+4.4 
-3.7 


1.11 


11.00 : 20.00 


3.60x10-3 


12 


+5.3 
-6.1 


1.11 



Table 1: Differential cross section for D production in pff , in the kinematic 
range 5 < < 1000 GeV'^ , 0.02 <y < 0.7, 1.5 < p^* < 20 GeV, \r]^'\ < 1.5. The 
columns list the bin range, the bin-averaged differential cross section, the statistical, 
^stat, o^nd systematic, Ssyst, uncertainties, and the QED correction factors, Cr- The 
overall normalization uncertainties from the luminosity (±1.9%) and branching 
ratio of the D*^ decay channel (±1.5%) are not included. 
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dcr 
dr)D* 

(nb) 


^stat 
(%) 


^syst 
(%) 


Cf 


-1.50 : -1.25 


1.48 


7.5 


+6.8 
-6.7 


1.06 


-1.25 : -1.00 


1.66 


5.4 


+5.6 
-5.3 


1.05 


-1.00 : -0.75 


1.61 


4.9 


+6.1 
-4.4 


1.05 


-0.75 : -0.50 


1.85 


4.2 


+4.6 
-3.8 


1.03 


-0.50 : -0.25 


1.94 


4.2 


+4.3 
-3.5 


1.03 


-0.25 : 0.00 


2.02 


4.0 


+4.3 
-3.7 


1.04 


0.00 : 0.25 


1.90 


4.4 


+4.2 
-3.4 


1.04 


0.25 : 0.50 


1.97 


4.4 


+4.3 
-3.3 


1.05 


0.50 : 0.75 


1.96 


4.7 


+4.5 
-3.6 


1.03 


0.75 : 1.00 


2.02 


4.9 


+4.8 
-4.2 


1.02 


1.00 : 1.25 


2.00 


5.8 


+5.3 
-5.1 


1.01 


1.25 : 1.50 


1.84 


7.7 


+7.4 
-5.6 


1.01 



Table 2: Differential cross section of D*^ production in rj^* . See TableUl for 
other details. 





dcr 
dzO* 

(nb) 


^stat 
(%) 


(^syst 
(%) 




0.000 : 0.100 


3.02 


12 


+8.6 
-7.1 


1.00 


0.100 : 0.200 


6.83 


6.1 


+6.1 
-5.0 


1.01 


0.200 : 0.325 


8.18 


3.5 


+5.5 
-4.9 


1.02 


0.325 : 0.450 


9.20 


2.5 


+4.6 
-3.8 


1.03 


0.450 : 0.575 


9.14 


2.3 


+4.6 
-4.0 


1.05 


0.575 : 0.800 


5.12 


2.4 


+6.5 
-5.1 


1.07 


0.800 : 1.000 


0.063 


9.1 


+9.9 
-8.5 


1.07 



Table 3: Differential cross section of D*"^ production in z^* . See Table U\ for 
other details. 
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( GeV2) 


da 

(nb/ GeV2) 


^stat 
(%) 


•^syst 
(%) 




5.0 : 8.0 


0.499 


3.9 


+6.7 
-6.1 


1.03 


8.0 : 10.0 


0.307 


4.3 


+6.0 
-5.2 


1.03 


10.0 : 13.0 


0.222 


4.0 


+4.9 
-4.1 


1.02 


13.0 : 19.0 


0.125 


3.5 


+5.6 
-5.0 


1.03 


19.0 : 27.5 


0.752x10-1 


3.7 


+4.9 
-4.0 


1.04 


27.5 : 40.0 


0.415x10-1 


3.9 


+4.8 
-3.8 


1.04 


40.0 : 60.0 


0.169x10-1 


4.7 


+5.6 
-5.6 


1.05 


60.0 : 100.0 


0.747x10-2 


5.0 


+7.1 

-5.1 


1.06 


100.0 : 200.0 


0.171x10-2 


7.8 


+6.6 
-4.4 


1.07 


200.0 : 1000.0 


0.140x10-3 


13 


+6.1 
-5.2 


1.14 



Table 4: Differential cross section of D*"^ production in . See TableU\for other 
details. 



y 


da 
dy 

(nb) 


"^stat 
(%) 


'^syst 
(%) 




0.02 : 0.05 


1.20 xlOi 


7.9 


+16 
-12 


1.07 


0.05 : 0.09 


2.07 xlOi 


3.4 


+6.7 
-6.5 


1.05 


0.09 : 0.13 


1.79 xlOi 


3.4 


+4.5 
-4.0 


1.04 


0.13 : 0.18 


1.37 xlOi 


3.6 


+4.6 
-4.8 


1.04 


0.18 : 0.26 


1.13 xlOi 


3.3 


+4.8 
-3.7 


1.04 


0.26 : 0.36 


8.03 


3.7 


+4.8 
-4.0 


1.03 


0.36 : 0.50 


5.09 


4.2 


+5.2 
-4.5 


1.02 


0.50 : 0.70 


2.90 


6.0 


+9.3 
-7.1 


1.01 



Table 5: Differential cross section of D*^ production in y. See Table\J\for other 
details. 
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X 


da 
dx 

(nb) 


<^stat 
(%) 


(%) 




(0.8 : 4.0) xlO"^ 


0.475 xlO^ 


3.5 


+6.0 
-5.3 


1.06 


(0.4 : 1.6) xlO-3 


0.198 xlO^ 


2.1 


+4.8 
-3.9 


1.03 


(1.6 : 5.0) xlO-3 


0.357 xlO^ 


2.6 


+4.9 
-3.9 


1.02 


(0.5 : 1.0) xlO-2 


0.553 xlO^ 


5.7 


+6.3 

-5.1 


0.99 


(0.1 : 1.0) xlO-i 


0.159 xlO^ 


10.7 


+9.2 
-8.4 


1.08 



Table 6: Differential cross section of D*^ production in x. See TableU\for other 
details. 
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y 


^ VIS 


'-'stat 


'-'syst 


beauty 


Cr 


( GeV^l 

I V.' V ' V J 


















020 
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1 20 

X ^ w 


23 


+19 

-20 




W . W 
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X . Wt: 




050 


• OQO 
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1 

X w 


+11 
-11 


1 5 

X . O 


1 04 

X . Wt: 


5 : 9 


OQO 


• 1 60 


421 


6 

w . w 


+6.8 
-7.0 


5 2 

W . ^ 


1 04 

X . Wt: 




0.160 


: 0.320 


550 


5.3 


+6.5 
-5.8 


11.0 


1.03 




320 


• 700 


456 


6 8 


+6.3 
-5.5 


18.2 


1.02 




020 


■ 050 


1 08 

± W(J 


14 

X^ 


+17 
-12 


1 

W . X 


1 05 

X . w w 




050 


• OQO 


1 78 


6 5 

W . W 


+7.0 
-6.0 


1 2 

X . ^ 


1 04 

X . Wt: 




OQO 


• 160 


220 


5 8 

W . (J 


+4.7 
-4.6 


2 9 


1 03 

X . w w 




0.160 


: 0.320 


352 


5.1 


+4.5 
-3.7 


8.1 


1.02 




0.320 


: 0.700 


307 


7.2 


+6.6 
-5.0 


12.5 


1.00 




020 


• 050 


65 1 

WW . ± 


1 5 

X w 


+ 16 
-12 


2 


1 07 

X . W 1 




050 


• OQO 


1 60 

± WW 


6 4 


+6.2 
-7.2 


1 2 

X . ^ 


1 04 

X . Wt: 


14 : 23 


OQO 


• 1 60 


205 


5 6 

W .U 


+4.7 
-4.7 


3 1 

O. X 


1 03 

X .WO 




0.160 


: 0.320 


267 


5.9 


+4.9 
-4.4 


9.0 


1.03 




320 


• 700 

\J • 1 \J\J 


250 

iii WW 


7 4 


+5.7 
-6.7 


13 5 

X w . w 


1 01 

X . W X 




020 


■ 050 


37 1 


29 


+ 18 
-18 


1 

W . X 


1 08 

X . Ww 




050 


• OQO 


1 34 


7 

1 .W 


+7.5 
-7.8 


9 

W . u 


1 06 

X . w w 


23 ■ 45 


OQO 


• 1 60 


1 Q6 

X <JW 


5 3 

W . W 


+4.4 
-4.3 


3 6 

W . W 


1 05 

X . w w 




0.160 


: 0.320 


275 


5.1 


+4.1 
-3.4 


10.2 


1.03 




0.320 


: 0.700 


284 


6.1 


+6.4 
-4.5 


14.7 


1.02 




0.020 


: 0.050 


14.2 


38 


+35 
-18 


0.0 


1.25 




0.050 


: 0.090 


72.1 


9.6 


+8.0 
-7.2 


1.2 


1.07 


45 : 100 


0.090 


: 0.160 


87.0 


8.4 


+4.9 
-4.6 


3.9 


1.04 




0.160 


: 0.320 


182 


5.7 


+5.3 
-3.9 


9.4 


1.04 




0.320 


: 0.700 


175 


7.6 


+6.6 
—5.6 


14.0 


1.02 




0.020 


: 0.350 


80.2 


11 


+7.6 
-4.2 


5.8 


1.1 


100 : 158 


0.350 


: 0.700 


45.1 


16 


+7.6 
-7.8 


5.0 


0.99 




0.020 


: 0.300 


59.8 


14 


+4.8 
-6.3 


3.5 


1.16 


158 : 251 


0.300 


: 0.700 


37.3 


17 


+6.6 
-4.9 


4.3 


1.04 




0.020 


: 0.275 


28.4 


24 


+8.2 

-10 


2.4 


1.26 


251 : 1000 


0.275 


: 0.700 


46.7 


21 


+8.7 
-5.1 


6.9 


1.07 



Table 7: Visible cross sections, a^is, for D*"^ production in bins of andy. The 
second but last column reports the estimated contribution from beauty decays, based 
on the Rapgap beauty MC rescaled to ZEUS data. See Table\J\for other details. 
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(GgV^) 




y 


<5l 


<52 




<53 


S4 




<55 


<56 


5, 






■Sll 


<5l2 


<5l3 




<^15 


<5l6 


<5l7 










(%) 


(%) 


(%) 


(%) 


(%) 


(%) 


(%) 


(%) 


(%) 


(%) 


(%) 


(%) 


(%) 


(%) 


(%) 
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+ 11 


+ 1 
+ 


8 
6 


+ 7 

— 2 
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-7.0 


-3.4 


+ 2 

— 2 


5 
4 


+ 2.5 


+ 8 

— 8 


8 
8 


+ 

— 






+ 
+ 






-4 
+4 


2 
2 


+ 1.8 


+ 

— 


1 
1 


+2.2 






0.050 


: 0.090 


+3.5 


+ 6 
— 4 


1 
8 


+ 2 
— 3 


3 



+6 
— 6 


6 
6 


+0.7 


-1.9 


+ 2 
— 2 


4 
2 


+ 2.1 


+ 3 
— 3 


4 
4 


-0 
+0 


1 

1 


+ 
+ 






-1 
+ 2 


9 
3 


+3.0 


-0 
+0 


7 
8 


-1.6 


5 


9 


0.090 


: 0.160 


+ 1.5 


+ 3 
-4 


3 



+ 2 
-1 


6 
5 


+ 1 
-1 


1 
1 


+ 1.0 
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+1 
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5 
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7 
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2 
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6 
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+ 


3 


+ 


7 
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7 
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6 
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-2 
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+ 1 
-1 


7 
7 
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2 
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1 


+0.8 


-1 

+1 


6 
7 
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0.160 
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+ 


2 

5 
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— 


3 
7 
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+0 
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5 


+0.3 


+0.5 


+ 2 


1 
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-1 


5 

5 


-0 
+0 


1 


+0 
— 


1 


+0 
— 


6 
7 
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-1 
+ 1 


8 
9 
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0.320 


: 0.700 
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+ 3.8 


+1 

-0 


1 
7 


+0 


1 


-0 
+0 


4 


+0.4 


-0.2 


+ 1 


9 
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+ 1 
-1 
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9 


+0 


2 


+0 


8 


+0 


4 
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-1 
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8 
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1 
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9 

8 


-1.7 


14 


23 
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+0 


4 
8 
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-0 


2 
3 
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+0 


3 
3 
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-1.6 
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-1 



9 
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+ 1 
-1 


6 
6 
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+ 
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1 
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-0 
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4 
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2 
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0.160 
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4 
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7 
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: 0.320 
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-0.1 
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+0.5 
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5 
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-0 
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7 


+0 
-0 






-0 
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+3.0 
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+ 2 
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1 
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0.320 
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-2.8 
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-0 
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3 


+0 
-0 
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-2.4 
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5 
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-2.2 


+4.4 


+ 1 
-1 


3 
3 


+ 1.4 
-1.4 


+ 1 
-1 


8 
8 


+ 0.9 

-1.1 


+0.0 
-0.0 


-0.2 
+ 0.2 


+ 1.3 


-1 
+ 1 


3 
2 


-0.2 
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0.350 


: 0.700 


-5.3 
+0.8 


-2 
-0 


5 
9 


+ 0.0 
-0.0 


-0.0 
+ 0.0 


+ 1.8 


+4.4 


+ 1 
-1 


4 
4 


+ 1.5 
-1.5 


+ 2 
-2 


9 
9 
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-3 
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2 


+0.5 
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+ 
-0 


1 
4 


-0.2 


-1 
+ 1 


3 
1 


+2.9 






0.020 


: 0.300 


+0.3 
-0.8 


-4.3 
+ 2.0 


+0.0 
-0.0 


-0.0 
+0.0 


-1.9 


-0.3 


+ 1 
-1 


2 
2 


+ 1.3 
-1.3 


+ 3.0 
-3.0 
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8 
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-0.6 
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-0 
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8 


+0.1 
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-2.6 
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-0.0 


-0.0 
+0.0 


+ 1.8 


+ 1.7 


+ 1 
-1 


4 
3 


+ 1.4 
-1.4 


+ 3.0 
-3.0 


+0 
-0 


1 
7 


+0.5 
-1.0 


+0.3 
-0.5 


+ 1.1 


-1.0 
+ 0.9 


+3.0 






0.020 


: 0.275 


-0.6 
-2.2 


-5.8 
+ 5.9 


+0.0 
-0.0 


-0.0 
+0.0 


-4.6 


-0.0 


+ 1 
-1 


2 
2 


+ 1.3 
-1.3 


+4.6 
-4.6 


+0 

-1 


8 
3 


+ 0.0 
+ 0.0 


-1.0 
+0.9 


-3.3 


+ 
-0 


1 
1 


-0.3 


251 : 


1000 


0.275 


: 0.700 


-0.8 
+4.3 


-1 
+ 3 


7 

8 


+0.0 
-0.0 


-0.0 
+0.0 


+2.5 


+3.5 


+ 1 

-1 


4 
3 


+ 1.4 
-1.4 


+ 3.6 
-3.6 


-0 
-0 


1 

1 


+0.7 

-1.4 


-0.5 
+0.4 


-1.4 


-1.0 
+ 0.9 


+0.9 



Table 8: Individual systematical uncertainties as defined in Section^ for the 
double- differential cross sections in bins of and y. The uncertainty 6s and 
6io are not reported as Sg is constant (+2%) and 5iq was found to be negligible. 
The overall normalisation uncertanties 6is = ±1.9% and 6ig = ±1.5% are also not 
listed. 
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5h 
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(%) 


(%) 


(%) 


(%) 


(%) 
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+0.2 
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+7.8 
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+0.6 
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7 
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+2.0 
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+5.0 


+0.1 
+0.8 
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+3.3 
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+0.2 


-0.9 
+ 1.2 
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-2.5 
+2.2 


±1.0 
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+5.7 


+4.7 
—6.3 


-1.6 
+2.5 


-2.4 
+5.9 


-4.0 
+4.2 


±2.2 
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+9.5 
—6.2 


-6.5 
+ 15 


+ 1.4 
+0.0 


-3.6 
+8.0 


+ 1.6 
+0.1 


±0.0 
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+0.1 
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+0.0 
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+ 13 
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+0.7 
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±0.1 
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+0.3 
—0.9 


-5.7 
+7.0 


+0.4 
—0.6 


-3.2 
+3.7 


-0.3 
— 0.6 


±0.4 
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+6.1 


+0.6 
+0.1 


-2.4 
+4.8 


-0.5 
+0.6 


±0.8 




0.00080 


-1.5 
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+3.2 


+0.3 
+0.3 


-1.9 

+4.4 


-0.9 
+0.7 


±1.8 
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+2.7 
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—3.7 


-1.0 
+ 1.0 


-2.5 
+4.5 
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±3.0 
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+0.6 
—0.5 
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+8.4 
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+3.9 
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1 A C 
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-|-U. 1 
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-l-o.y 
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-0.4 
— 1.6 
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+0.4 


-1.8 
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+ 1.6 
+0.4 


±0.0 
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1 7 


-4.8 

-|-D.U 


-0.3 
n 7 


-1.9 


-0.1 

— U.D 


±0.9 


60 
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-0.2 
+0.8 


-3.9 
+5.2 


+0.1 
-0.0 


-1.4 
+2.7 
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+0.3 


±2.4 
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-0.9 
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-0.2 
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+0.0 


±2.9 
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±4.8 




0.01000 
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+5.3 


+0.4 
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-1.5 
+2.3 


+0.0 
+0.3 


±4.2 
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-0.8 
+ 1.3 


-2.0 
+2.3 


+0.4 
-0.5 


-1.3 
+ 1.9 


-1.0 
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±7.2 
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-0.1 
-0.1 


-3.7 
+3.8 


+0.4 
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+ 1.4 
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±4.1 


200 


0.00500 
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-0.6 
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-0.0 


-0.7 
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±5.1 


350 


0.01000 


-0.2 
+ 1.3 


-2.8 
+3.7 


+0.0 
+0.3 


-0.6 
+0.9 


+0.0 
+0.1 


±10.6 



Table 9: Breakdown of the theoretical uncertainty on (T^g^(x, Q^), showing the 
uncertainty from the variation of charm mass (Smc), of the renormalisation and 
factorisation scales (S^), of as (^as)> of the fragmentation function (Saj^), of the 
transverse fragmentation (Skj.), and of the expected beauty component (6]^). The 
upper (lower) value gives the effect of a positive (negative) variation of the para- 
meter. 
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X 
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A 
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oooso 
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1 24 
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7 
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00030 


0.187 


5.5 
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14 
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12 
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6 
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00050 
00022 
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7 8 
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-3.8 
+6.9 
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+6.0 
-4.9 
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-7.2 


0.492 
332 


IS 


00450 
00250 
001 35 


081 
1 68 
1 99 


15 

X O 

6 5 
5 8 


+ 16 
-12 
+6.2 
-7.2 
+4.7 
-4.8 


+17 
-9.5 
+8.0 
-6.7 
+7.9 
-5.1 


0.262 
456 
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00080 


0.216 


6.3 


+5.1 
-4.6 


+5.9 
-5.0 


0.544 




00035 


324 


8 3 


+6.1 
-7.1 


+6.5 
-6.8 


370 




0.00800 


0.068 


29 


+ 18 
-19 


+15 
-11 


0.236 




00550 


1 59 


7 1 

1 . X 


+7.5 
-7.9 


+9.1 
-6.2 


478 


32 


00240 


234 


5 5 


+4.5 
-4.4 


+7.6 
-4.1 


560 




00140 


0.266 


5.5 


+4.3 
-3.6 


+6.5 
-4.4 


0.594 




00080 


0.390 


6.8 


+6.8 
-4.7 


+5.7 
-5.5 


0.430 




01500 


0.068 


38 


+35 
-18 


+15 
-8.6 


0.166 


60 


00800 
00500 


0.173 
0.161 


9.9 
9.2 


+8.1 
-7.3 
+5.2 
-4.9 


+6.6 
-5.6 
+6.4 
-4.8 


0.465 
0.570 




00320 

vy • vy vy ^— ' w 


0.258 


6.4 


+5.6 
-4.1 


+6.6 
-5.1 


0.624 




0.00140 


0.327 


9.0 


+7.2 
—6.2 


+6.0 
—6.0 


0.516 




0.01000 


0.130 


12 


+8.2 
-4.6 


+7.2 
-6.4 


0.491 


120 


0.00200 


0.287 


21 


+8.7 
-8.9 


+7.9 
-7.7 


0.583 




0.01300 


0.176 


17 


+5.2 
-6.8 


+5.8 
-5.6 


0.465 


200 


0.00500 


0.239 


22 


+7.6 
-5.6 


+8.6 
-8.8 


0.624 




0.02500 


0.103 


29 


+9.0 

-11 


+6.2 
-6.4 


0.433 


350 


0.01000 


0.193 


29 


+ 11 
-6.2 


+11 
-11 


0.636 



Table 10: The reduced cross-section (T^g^(a:, Q^) with statistical, systematic and 
theoretical uncertainties. The last column shows the kinematical acceptance. 
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Figure 1: Distribution of M{Kir) for D*^ candidates with 143.2 < AM < 
147.7 MeV (filled circles) and for wrong-sign combinations (empty circles). The 
signal region is marked as a shaded area. 
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Figure 2: Distribution of the mass difference, AM = M^KtttIs) — M{Ktt), 
for the D*^ candidates with 1.80 < M{Kit) < 1.92 GeV (filled circles) and for 
wrong-sign combinations (empty circles). The background fit described in the text 
is shown as a dashed (continuous) line for correct-sign (wrong-sign) combinations. 
The D*^ signal region is marked as a shaded area. 
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Figure 3: Number of reconstructed D*^ (filled circles), divided by bin size, as a 
function of , t]^* , Q^, y and . Data are compared to a MC mixture contain- 
ing non-diffractive and diffractive charm production in DIS, beauty production, and 
charm photoproduction. The sum of the MC samples is normalised to the number 
of D*^ in the data. 
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Figure 4: Differential D*"^ cross sections as a function of (a) , (h) rj^* and 
(c) z^* (filled circles). The error bars show the statistical and systematic uncer- 
tainties added in quadrature, the inner bars show the statistical uncertainties alone. 
Also shown are NLO QCD predictions calculated using HvQDis (dashed line and 
shaded area for the uncertainties) and Rapgap MC prediction for charm creation 
via boson-gluon fusion (long-dashed line). The contribution from b-quark decays, 
calculated with the Rapgap MC (continuous line), is included in the predictions. 
The MC cross sections for charm (beauty) are scaled by 1.1 (l.Q) as described in 
the text. 
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Figure 5: Differential D cross sections as a function of (a) Q , (h) y and (c) 
X. Other details as in Fig. [^j 
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Figure 6: D ouble- differential D*^ cross sections as a function of and y for 
5 < < 100 GeV"^ (filled circles). The measurements from the HI collaboration 
(empty squares) are also shown fl^ . Other details as in Fig. [7} 



29 



ZEUS 



CM 
> 

O 
a. 

CM 

o 

■a 
>« 
;o 

CM 



CM- 0-3 

« 0.25 
CD 

0.2 
c^ 0.15 

o 

■o 0.1 

CM 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 

y 



' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I 

Q2=251-1000 GeV^ 



I ■ ■ ■ ■ I ■ ■ ■ ■ I ■ ■ ■ ■ I ■ ■ ■ ■ I ■ ■ ■ ■ I 



CM 
> 

O 

!a 
a. 

CM 

o 

>. 

:a 

CM 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 

y 

ep ^ e D* X 

ZEUS D* 363 pb ^ 
A H1 D* (high Q^) 
[ HVQDIS + RAPGAPbx1.6 

RAPGAPBGFcxl.1 +bx1.6 

RAPGAPbx1.6 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 



Figure 7: D ouble- differential D*^ cross sections as a function of and y for 
100 < < 1000 GeV"^ (filled circles). The measurements from the HI collabora- 
tion (empty triangles) are also shown [T^ . Other details as in Fig. \^ 



30 



ZEUS 



0) 

-§ 1.8 

^1.4 
1.2 
1 

0.8 
0.6 
0.4 
0.2 




T 1 1 1 — I — I — n~ 



T 1 1 1 — I — I — r~r 



• ZEUS D* 363 pb 



-1 



I I 



10 



J I I I I I L_l_ 



10^ 



J I I I I I L_L 



10^ 
(GeV^) 



Figure 8: Ratio of e p to e~^p visible D*^ cross sections as a function of Q'^ . 
Only statistical uncertainties are shown. Bin boundaries are as in Table^ 
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Figure 9: Reduced charm cross sections from D*"^ (filled circles) compared to 
the ZEUS measurement (empty squares) and the combination of previous 
HERA results (empty circles). The outer error bars include experimental and 
theoretical uncertainties added in quadrature. The inner error bars in the ZEUS 
D* and measurements show the experimental uncertainties. The inner error 
bars of the combined HERA data represent the uncorrected part of the uncertainty. 
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Figure 10: Reduced charm cross sections (filled circles) compared to a GM-VFNS 
calculation based on HERAPDF1.5 parton densities. The inner error bars show the 
experimental uncertainties and the outer error bars show the experimental and the- 
oretical uncertainties added in quadrature. The outer bands on the HERAPDF1.5 
predicition show the total uncertainty while the inner bands correspond to the sum 
in quadrature of all uncertainties excluding the charm-mass variation. 
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